Three-dimensional stability of Burgers vortices 



Thierry G allay 
Institut Fourier 
Universite de Grenoble I 
BP 74 

38402 Saint-Martin-d'Heres, France 
Thierry . GallayOuj f -grenoble . f r 



1-1 Rokkodai, Nada-ku 
Kobe 657-8501, Japan 
yasunor iOmath . kobe-u .ac.jp 



Yasunori Maekawa 
Faculty of Science 
Kobe University 



February 9, 2010 



Abstract 



Burgers vortices are explicit stationary solutions of the Navier-Stokes equations which 
are often used to describe the vortex tubes observed in numerical simulations of three- 
dimensional turbulence. In this model, the velocity field is a two-dimensional perturbation of 
a linear straining flow with axial symmetry. The only free parameter is the Reynolds number 
Re = r/i/, where T is the total circulation of the vortex and v is the kinematic viscosity. The 
purpose of this paper is to show that Burgers vortex is asymptotically stable with respect 
to general three-dimensional perturbations, for all values of the Reynolds number. This 
definitive result subsumes earlier studies by various authors, which were either restricted to 
small Reynolds numbers or to two-dimensional perturbations. Our proof relies on the crucial 
observation that the linearized operator at Burgers vortex has a simple and very specific 
dependence upon the axial variable. This allows to reduce the full linearized equations 
to a vectorial two-dimensional problem, which can be treated using an extension of the 
techniques developped in earlier works. Although Burgers vortices are found to be stable 
for all Reynolds numbers, the proof indicates that perturbations may undergo an important 
transient amplification if Re is large, a phenomenon that was indeed observed in numerical 
simulations. 

1 Introduction 

The axisymmetric Burgers vortex is an explicit solution of the three-dimensional Navier-Stokes 
equations which provides a simple and widely used model for the vortex tubes or filaments 
that are observed in turbulent flows [UlSOj. Despite obvious limitations, due to oversimplified 
assumptions, this model describes in a correct way the fundamental mechanisms which are 
responsible for the persistence of coherent structures in three-dimensional turbulence, namely 
the balance between vorticity amplification due to stretching and vorticity dissipation due to 
viscosity. If one believes that vortex tubes play a significant role in the dynamics of turbulent 
flows, it is an important issue to determine their stability with respect to perturbations in the 
largest possible class. So far, this problem has been studied only for the axisymmetric Burgers 
vortex and for a closely related family of asymmetric vortices |27l [21] . 

As was shown by Leibovich and Holmes [19J, one cannot hope to prove energetic stability 
of the Burgers vortex even if the circulation Reynolds number is very small. To tackle the 
stability problem, it is therefore necessary to have a closer look at the spectrum of the linearized 



operator. This is a relatively easy task if we restrict ourselves to two-dimensional perturbations. 
Assuming that the vortex tube is aligned with the vertical axis, this means that the perturbed 
velocity field lies in the horizontal plane and does not depend on the vertical variable. Under 
such conditions, the Burgers vortex is known to be stable for any value of the Reynolds number. 
This result was first established by Giga and Kambe [15] for Re <C 1 and then by Gallay and 
Wayne [H] in the general case. Moreover, a lot is known about the spectrum of the linearized 
operator, which turns out to be purely discrete in a neighborhood of the origin in the complex 
plane. Using perturbative expansions, Robinson and Saffman [27j showed that all linear modes 
are exponentially damped for small Reynolds numbers. This property was then numerically 
verified by Prochazka and Pullin [25] for Re < 10^, and finally rigorously established in |llj . 

The situation is much more complicated if we allow for arbitrary three-dimensional pertur- 
bations. In that case, it was shown by Rossi and Le Dizes [28j that the linearized operator does 
not have any eigenfunction with nontrivial dependence in the vertical variable. While this result 
precludes the existence of unstable eigenvalues, it also implies that stability cannot be deduced 
from such a simple analysis, and that continuous spectrum necessarily plays an important role. 
Unfortunately, the vertical dependence of the perturbed solutions is not easy to determine, as 
can be seen from the note [3] where a few attempts are made in that direction. The only rigorous 
result so far is due to Gallay and Wayne [12], who proved that the Burgers vortex is asymptoti- 
cally stable with respect to three-dimensional perturbations in a fairly large class provided that 
the Reynolds number is sufficiently small. For larger Reynolds numbers, up to Re = 5000, an 
important numerical work by Schmid and Rossi [29] indicates that all modes are exponentially 
damped by the linearized evolution, although significant short-time amplification can occur. 

In this paper, we prove that the axisymmetric Burgers vortex is asymptotically stable with 
respect to three-dimensional perturbations for arbitrary values of the Reynolds number. As in 
[12] . we assume that the perturbations are nicely localized in the horizontal variables, but we 
do not impose any decay with respect to the vertical variable. Our approach is based on the 
fact that the linearized operator has a very simple dependence upon the vertical variable: the 
only term involving X3 is the dilation operator x^dxg, which originates from the background 
straining field. This crucial property was already exploited in [281 [3l [29]. but we shall show that 
it allows to reduce the three-dimensional stability problem to a two-dimensional one, which can 
then be treated using an extension of the techniques developped in [11]. Although the spectrum 
of the linearized operator remains stable for all Reynolds numbers, the estimates we have on 
the associated semigroup deteriorate as Re increases, in full agreement with the amplification 
phenomena observed in [29] . 

We now formulate our results in a more precise way. We start from the three-dimensional 
incompressible Navier-Stokes equations: 

dtV + (V,V)V = uAV --VP , V-V = 0, (1.1) 

P 

where V = V{x, t) E M"^ denotes the velocity field, P = P{x, t) G M is the pressure field, and 
X = {xi, X2, X3)~^ G is the space variable. The parameters in (|1.1|) are the kinematic viscosity 
u > and the density p > 0. To obtain tubular vortices, we assume that the velocity V can be 
decomposed as follows: 

V{x,t) = V'{x) + U{x,t) , (1.2) 
where is an axisymmetric straining fiow given by the explicit formula 

/-xi\ /-i 0\ 

V'{x) = ^ -X2 = jMx , where M = . (1.3) 

V2X3/ V 1/ 
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Here 7 > is a parameter which measures the intensity of the strain. Note that V -V^ = 0, and 
that is a stationary solution of (II. ip with the associated pressure = — |^. Our goal 
is to study the evolution of the perturbed velocity field U{x,t). 

To simplify the notations, we shall assume henceforth that 7 = = p = 1. This can be 
achieved without loss of generality by replacing the variables x, t and the functions V, P with 
the dimensionless quantities 

7\i/2 , . V ~ P 

^ X . t = lt, V = —— , P 



(71/) 1/2 ' P71/ 

For further convenience, instead of considering the evolution of V or U, we prefer working with 
the vorticity field Q = V xV = V xU . Taking the curl of and using (|1.2|) . (|1.3|) . we obtain 
for 17 the evolution equation 

dt^ + {u,v)n - {n,V)u = LQ , v-n = o, (i.4) 

where L is the differential operator defined by 

in = An - {Mx, + Am . (1.5) 

Under mild assumptions that will be specified below, the velocity field U can be recovered 
from the vorticity il. via the three-dimensional Biot-Savart law 

= _i_ / '--""'"fa' d, =: iK,o,n)W) . (1.6) 

In what follows we shall often encounter the particular situation where the velocity U is two- 
dimensional and horizontal, namely U{x) = {Ui{xh),U2{xh),0)^ where Xh = {xi,X2)^ € M^. 
In that case the vorticity satisfies Q{x) = {0,0,Q,^(xh))~^ , and the relation (|1.6|) reduces to the 
two-dimensional Biot-Savart law 



Uhixh) = ^ [ ^-^^-^^nMdVh =: {K2D^^z){xh) , (1.7) 
2-K J^2 \xh-yhr 




where Uh = {Ui,U2)^ and = {—X2,xi)^. 

We can now introduce the Burgers vortices^ which are explicit stationary solutions of ([T7 
of the form Vt = aG, where a € M is a parameter. The vortex profile is given by 

G{x) = I I , where ^(x,,) = ^ e^l^'^l'/^ . (1.8) 

The associated velocity field U = all^ can be obtained from the Biot-Savart law (|1.7p and has 
the following form 

U^(x) = u^dxhl"^) i xi] , where u9{r) = -i- (^1 - e"^/^) . (1.9) 

If J7 = aG, it is easy to verify that a = ^six^) dxh- This means that the parameter a S M 
represents the total circulation of the Burgers vortex aG. In the physical literature, the quantity 
\a\ is often referred to as the (circulation) Reynolds number. 
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The aim of this paper is to study the asymptotic stabiHty of the Burgers vortices. We thus 
consider solutions of (II. 4p of the form = aG + uo^ U = aU'^ + n, and obtain the following 
evolution equation for the perturbation: 

+ (n, V)a; - (w, V)u = (L - qA)w , V • w = , (1.10) 

where A is the integro-differential operator defined by 

kuj = (C/^, V)a; - (cj, V)U^ + (n, V)G - (G, V)n . (1.11) 

Here and in the sequel, it is always understood that u = K^d * w. 

An important issue is now to fix an appropriate function space for the admissible pertur- 
bations. Since the Burgers vortex itself is essentially a two-dimensional flow, it is natural to 
choose a functional setting which allows for perturbations in the same class, but we also want 
to consider more general ones. Following |I2] . we thus assume that the perturbations are nicely 
localized in the horizontal variables, but merely bounded in the vertical direction. As we shall 
see below, this choice is more or less imposed by the particular form of the linear operator (jl.5p . 

To specify the horizontal decay of the admissible perturbations, we first introduce two- 
dimensional spaces. Given m € [0, oo], we denote by pm ■ [0, oo) — > [l,oo) the weight function 
defined by 

1 if m = , 

Pmir) = + if 0<m<(X), (1.12) 

^ e^/^ if m = oo . 

We introduce the weighted space 

L\m) = { f £ L^R^) [ \f{xh)\^pm{\xh\^)dxh<(x>} , (1.13) 

which is a Hilbert space with a natural inner product. Using Holder's inequality, it is easy to 
verify that Lp'{m) ^ L^(E?) if m > 1. In that case, we also define the closed subspace 

Ll{m) = {/ G L\m) [ fixh)dxh = o| . (1.14) 

Next, we define the three-dimensional space X{m) as the set of all : M'^ ^ M for which 
the map Xh i-> (j){xh,xs) belongs to L?'{m) for any x^ G M, and is a bounded and continuous 
function of X3. In other words, we set 

X{m) = BC{R;L^{m)) , Xo(m) = 5C(M;Lg(m)) , (1.15) 

where "i?C(R;y)" denotes the space of all bounded and continuous functions from M into Y. 
Both X{m) and Xo(m) are Banach spaces equipped with the norm 

ll'/'llx{m) = sup ||(/'(-,2;3)||l2(^) . (1.16) 

X3GR 

Our goal is to study the stability of the Burgers vortex = aG with respect to perturbations 
oj € X{mY . In fact, we can assume without loss of generality that uj belongs to the subspace 

X(m) = X{m) X X{m) x Xo(m) C X{mf , (1.17) 

which is invariant under the evolution defined by (|1.1U|) . This is a consequence of the following 
result, whose proof is postponed to Section [6?T] 
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Lemma 1.1 Fix m G (l,oo]. If Cj £ X(in)^ satisfies V • = in the sense of distributions, 
then there exists d G M such that 

/ Ld3{xh,X3)dxh = a , forallx3eR. (1-18) 

In view of Lemma ll. 11 if = aG + Co for some Cj E X{m)^ , we can write Q. = {a + a)G + w, 
where a is given by (jl.lSp and u = Co — aG. Then uj E X(m) by construction, and we are led 
back to the stabihty analysis of the Burgers vortex (q + a)G with respect to perturbations in 
X(m). 

In what follows we always consider the solutions uj{x,t) of (ll.lOp as X(m)-valued functions 
of time, and we often denote by uj{-,t) or simply uj(t) the map x i— )• uj{x,t). A minor drawback 
of our functional setting is that we cannot expect the solutions of (ll.lOp to be continuous in 
time in the strong topology of X(m). This is because the operator L defined in ()1.5p contains 
the dilation operator —x^dx^, see Section [2T] below. To restore continuity, it is thus necessary 
to equip X(m) with a weaker topology. Following [12], we denote by Xiodm) the space X{m) 
equipped with the topology defined by the family of seminorms 

\\(l>\\x„(m) = sup |l0(-,X3)||i2(^) , nGN. 

1 <n 

In analogy with (I1.17p . we set Xioc{rn) = Xioc{m) x Xioc{m) x Xo,/oc("i), where Xo^iodrn) is of 
course the space Xo(m) equipped with the topology of Xiodm). 

We are now able to formulate our main result: 

Theorem 1.2 Fix m G (2, oo] and a G M. Then there exist 5 = 5(q, m) > and G = G{a, m) > 
1 such that, for any loq G X(m) with V ■ ujq = and ||wo||x(m) ^ Eq. (ll.lOp has a unique 
solution UJ G L°°(M+ ;X(m)) PI C([0, oo) ;X;oc('n')) with initial data ojq. Moreover, 

\Ht)\\x(n.) < C||a;o||xM ' aUt>0 . (1.19) 

Theorem 11.21 shows that the Burgers vortex aG is asymptotically stable with respect to 
perturbations in X(m), for any value of the circulation a G M. If one prefers to consider 
perturbations in the larger space X(m)^, then our result means that the family {aG}aeR of all 
Burgers vortices is asymptotically stable with shift, because the perturbations may then modify 
the circulation of the underlying vortex. The key point in the proof is to show that the linearized 
operator L — aA has a uniform spectral gap for all a G M. This implies a uniform decay rate in 
time for the perturbations, as in (I1.19p . However, it should be emphasized that the constants 
G and 5 in Theorem 11.21 do depend on a, in such a way that C(a, m) — )• oo and 6{a, m) — )• as 
|a| oo. This is in full agreement with the amplification phenomena numerically observed in 

m- 

The proof of Theorem 11.21 gives a more detailed information on the solutions of (ll.lOp than 
what is summarized in (jl.l9p . First of all, we can prove stability for any m > 1, but the 
exponential factor e~*/^ in (I1.19P should then be replaced by e"''*, where r/ < (m — l)/2 if 
m < 2. Next, thanks to parabolic smoothing, we can obtain decay estimates not only for uj{t) 
but also for its spatial derivatives. Finally, due to the particular structure of the linear operator 
L — aA, it turns out that the horizontal part uj^ = {^1t^2)~^ of the vorticity vector has a faster 
decay than the vertical component as t — )• oo. Thus, a more complete (but less readable) 
version of our result is as follows: 
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Theorem 1.3 Fix m € (l,oo], a € M, and take € (1, |), rj G (0, ^] such that 2/z < m + 1 
and 2r] < m — 1. Then there exist 6 = 6{a,m) > and C = C{a,m, fi,T]) > 1 such that, for 
all initial data loq € X(m) with V • = and ||wo||x(m) ^ ^Q- (|1-10P has a unique solution 
uj G L°°(R+ ; X(m)) n C([0, oo) ; X;oc(r«)). Moreover, for all t > 0, 



C\\uJo\ 



MmxM^ < e-^' , (1.20) 



l|9fu.3(t)||xM<^;|^e-*, (1.21) 



where a{t) = 1 — e * and /3 G is any multi-index of length \f3\ = /3i + /32 + /^a < 1- 



The decay rates ()1.20p . ()1.2ip are optimal when /3 = 0, but it turns out that vertical deriva- 
tives such as dx^0Jh{t) or dx^ujz{t) have a faster decay as t — )• oo, see Sections H] and [S] for more 
details. In any case, we believe that the optimal rates are those provided by the linear stability 
analysis, as in Proposition 14. II below. 

The rest of this paper is devoted to the proof of Theorems 11.21 and 11.31 Before giving the 
details, we explain here the main ideas in an informal way. As was already mentioned, the main 
difficulty is to obtain good estimates on the solutions of the linearized equation 

dtoj = [L- ak)uj , V • cj = . (1.22) 

Once this is done, the nonlinear terms in (ll.lOp can be controlled using rather standard argu- 
ments, which are recalled in Section[5j To study p.22p . we use the fact that the operator L — qA 
depends on the vertical variable in a simple and very specific way. Indeed, it is easy to verify 
that [9^.3, L] = —dx^ and [92^3,A] = 0. This key observation, which already plays a crucial role 
in the previous works ^28i i3t i29j, implies the following identity: 

Qk^f,t{L-aK)^^ = e-^*e*(-^-"^)5^3a;o , (1-23) 

for all A; G N and all t > 0. If we take A; G N sufficiently large, depending on |a[, we can use (I1.23P 
to to show that d^^uj{t) decays exponentially as t ^ cxd if uj{t) is a solution of ()1.22p . Then, 
by an interpolation argument, we deduce that all expressions involving at least one vertical 
derivative play a negligible role in the long-time asymptotics, see Section H] for more details. 
This "smoothing effect" in the vertical direction is due to the stretching properties of the linear 
flow (fL2D . 

As a consequence of these remarks, we can restrict our attention to those solutions of (|1.22p 
which are independent of the vertical variable X3. We call this particular situation the vectorial 
2D problem, and we study it in Section [3l Note that the perturbations we consider here are 
two-dimensional in the sense that dx-j^u = dx^^u} = 0, but that all three components of u or a; 
are possibly nonzero. This is in contrast with the purely two-dimensional case considered in 

112) . where in addition U3 = uji = UJ2 = 0. Nevertheless, it is possible to show that the 
solutions of (I1.22P with dx^uj = converge exponentially to zero as t ^ cxd, and that the decay 
rate is uniform in a. Extending the techniques developped in |1H [12], this can be done using 
spectral estimates and a detailed study of the eigenvalue equation (L — aA)uj = Xu. It is then a 
rather straightforward task to complete the proof of Theorem 1 1 . 2 1 using the arguments presented 
above. 

Remark. The vortex tubes observed in numerical simulations are usually not axisymmetric: 
in general, they rather exhibit an elliptical core region. A simple model for such asymmetric 
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vortices is obtained by replacing the straining flow in (jl.3p with the nonsymmetric strain 
V^{x) = 'yM\x, where A € (0, 1) is an asymmetry parameter and 

l+A 

Ma = I 0^ -i^ I . (1.24) 

Asymmetric Burgers vortices are then stationary solutions to (II. 4p , where the operator L in the 
right-hand side is defined by (|1.5p with M replaced by Mx. Unlike in the symmetric case A = 0, 
no explicit formula is available and proving the existence of stationary solutions is already a 
nontrivial task, except perhaps in the perturbative regime where either the asymmetry parameter 
A or the circulation number a is very small. In view of these difficulties, asymmetric Burgers 
vortices were first studied using formal asymptotic expansions and numerical calculations, see 
6-g- |27 1 12 H [26]. The mathematical theory is more recent, and includes several existence results 
which cover now the whole range of parameters A E (0, 1) and a € M [12l[l3l|22l|23]. In addition, 
the stability with respect to two-dimensional perturbations is known to hold at least for small 
values of the asymmetry parameter [131122]. However, the only result so far on three-dimensional 
stability is restricted to the particular case where the circulation number a is sufficiently small, 
depending on A ■ 

Using Theorem 11.21 and a simple perturbation argument, it is easy to show that asymmetric 
Burgers vortices are stable with respect to three-dimensional pertubations in the space X(m), 
provided that the asymmetry parameter A is small enough depending on the circulation number 
a. This follows from the fact the the linearized operator at the symmetric Burgers vortex has 
a uniform spectral gap for all a G M, and that the asymmetric Burgers vortex is 0(A) close 
to the corresponding symmetric vortex in the topology of X(m), uniformly for all q € M |13] . 
Although this stability result is new and not covered by [12], it is certainly not optimal, and we 
prefer to postpone the study of the three-dimensional stability of asymmetric Burgers vortices 
to a future investigation. 



2 Preliminaries 

In this preliminary section we collect a few basic estimates which will be used throughout the 
proof of Theorems 11.21 and 11.31 They concern the semigroup generated by the linear operator 
(jl.Sp . and the Biot-Savart law (jl.6p relating the velocity field to the vorticity. Most of the 
results were already established in [HI Appendix A], and are reproduced here for the reader's 
convenience. 

As in [12] , we introduce the following generalization of the function spaces ()1.13p and ()1.15p . 
Given m € [0, oo] and p £ [1, oo), we define the weighted space 

LP{m) = {/ e LP(M2) |||/||^,(^) = \f{x,)rpU\xh\'r/'dx, < oo} , 

and the corresponding three-dimensional space 

XP{m) = BC{R;LP{m)) , H^IU^M = ^up ||0(-, X3) lU^M • 

If m > 2 — |, we also denote by L^{m) the subspace of all / € LP{m) such that f dx^ = 0. In 
analogy with (fTTTH . wesetXP(m) = XP (m) x XP (m) x (m) , where X^(m) = BC{R; Ll{m)). 
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2.1 The semigroup generated by L 

If we decompose the vorticity oj into its horizontal part uj^ = ("^ii ^^2)^ and its vertical component 
wa, it is clear from and (|1.5|) that the linear operator L has the following expression: 



Lu. = f^"^") = f ( V ^V^^^'''') > (2-1) 



where £/i is the two-dimensional Fokker-Planck operator 



= A, + :^.V. + 1 = E<+Ey^-.+1' (2-2) 



and £3 = — X3(9a;3 is a convection-diffusion operator in the vertical variable. 

As is shown in [TUl appendix A], the operator Ch is the generator of a strongly continuous 
semigroup in L?'{m) given by the explicit formula 

{e'^^ f){xh) = I e-^^^/(W/')dy/. , t>0, (2.3) 



A-Ka{t) 

where a{t) = 1 — e~*. Similarly, the operator £3 generates a semigroup of contractions in 
given by 

(e*^^/)(x3) = ^— — — / e 2.(2.) /(yg) , f > , (2.4) 
^2710(2^) Jm 

see [121 Appendix A]. Note that the semigroup e*^^ is not strongly continuous in the space 
-BC(M) equipped with the supremum norm. This is mainly due to the dilation factor e~* in 
(j2.4p . However, if we equip i3C(R) with the (weaker) topology of uniform convergence on 
compact sets, then the map 1 1— )■ e*^^/ is continuous for any / E i?C(M). This observation is the 
reason for introducing the space Xiodm) in Section [TJ 

Since the operators Ch and £3 act on different variables, it is easy to obtain the semigroup 
generated by L^^ = Ch + C3 by combining the formulas (|2.3p and ()2.4p . We find 

(e*^='(/))(x) = , ^ / e 5^^^^ e*^" </>(•, y3) (x/.)dy3 , t>0. (2.5) 

^y27^a{2t) Jr ^ ^ 

In [121 Proposition A. 6] , it is shown that this expression defines a uniformly bounded semigroup 
in X{m) for any m > 1, and that the map t 1— > e*^^ is strongly continous in the topology of 
Xioc{m). Moreover, the subspace Xo(m) is left invariant by e*^^ for any t > 0. Using these 
results and the relation (|2.ip . we conclude that the three-dimensional operator L generates a 
uniformly bounded semigroup in the space X(m), given by 



e 



*^a; = (e"3*/2e*-^3^i , e'^^/^e^^^ws , e^^'^s) , * > . (2.6) 



As is easily verified, if V • = 0, then V • e^^uj = for all t > 0. 

The asymptotic stability of the Burgers vortices relies heavily on the decay properties of the 
semigroup e*^ as t — )• 00. In the proof of Theorems 11.21 and 11.31 we also use the smoothing 
properties of the operator e*^ for t > 0, and in particular the fact that e*^ extends to a bounded 
operator from 'KP{m) into X^(m) for all p E [1, 2]. All the needed estimated are collected in the 
following statement. 
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Proposition 2.1 Let m G (l,oo], p G [1,2], and take rj G (0, ^] such that 21] < m — 1. For any 

P = /?2) /^s) G N^, there exists C > such that the following estimates hold: 

3 tL Ce-(i+^«)* 
a(t)p 2+ 2 

e ^W3||x(m) < i_i 1^ ll^3||xP(m.) > (2-8) 

a(t)p 2+^ 

for any oj G W(m) and all t > 0. Here a{t) = 1 — e^* and \ = /9i + /32 + /Ss. 

Proof. We first assume that m G (l,oo). If p G [1,2] and j3h = (/3i,/32) £ N^, it is proved in 
[TOj Appendix A] that 

l|5f^*''VllL2(^) < ^ II/IIlpM , t > , (2.9) 

a(t)f 2+ 2 

for all / G LP{m). If in addition / G LQ(m), we have the stronger estimate 

ll<^e*^''/llL2(^) < , t>0, (2.10) 

a(t)p 2+ 2 

where 77 > is as in Proposition 12. 11 On the other hand, using (j2.5p . we find by direct calculation 

llC'e^'^'/llL-w < — 3rll/lli-w ' * > • (2-11) 

a{t)- 

Here, as in (jl.23p . the stabilizing factor e~^^* comes from the dilation operator —x^dx^ which 



enters the definition of C3. Now, if we start from the representation ()2.5p and use the estimates 
()2.9p - ()2.1ip . we easily obtain ()2.7p . (|2.8p by a direct calculation, see [121 Proposition A. 6]. 

To complete the proof of Proposition 12. H it remains to show that (|2.9p , (|2.1Up still hold 
when m = 00. If t G (0,1), estimate (j2.9p is easily obtained by a direct calculation, based on 
the representation (j2.3p . Using this remark and the semigroup property of e*^'', we conclude 
that it is sufficient to establish (j2.9p . (|2.10p in the particular case where p = 2 and (3h = 0. 
This in turns follows easily from the spectral properties of the generator Ch- Indeed, it is 
well-known that Ch is a self-adjoint operator in L^(oo) with purely discrete spectrum cr{Ch) = 
{— I I /c = 0, 1, 2, . . . }. Moreover, the subspace Lq{oo) is precisely the orthogonal complement of 
the eigenspace corresponding to the zero eigenvalue, see for example [111 Lemma 4.7]. It follows 
that e*^'' is a semigroup of contractions in L^(oo), and that ||e*^''/||L2(oo) ^ e~*^^ll/llL2(oo) foi" 
alH > if / G Ll{oo). This proves ([221) and ([212]), with t/ = 1/2. □ 

2.2 Estimates for the velocity fields 

If the velocity u and the vorticity a; are related by the Biot-Savart law (jl.6p . we have < J(|tj[), 
where J is the Riesz potential defined by 

J{<P){x) = ^ [ ^^0(y)dy, xeR^. (2.12) 

Since lo will typically belong to the Banach space X(m), we need estimates on the Riesz potential 
J(0) for (j) G X(m). We start with a preliminary result: 
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Lemma 2.2 Let pi e [1,2), p2 £ [1,2], and assume that cp G X^i (0) n (0) . Ifqi,q2 G [l,oo] 
satisfy 

< qi < oo , P2 < q2 < , (2.13) 

2-pi 2-P2 

then J{(t)) = Ji{(t)) + J2{<P) with Ji{4>) € X'^'(O) for i = 1,2, and we have the following estimates 

II Ji (0)11x^1(0) < C'(pi,gi)||0||xPi(o) , (2.14) 
\\J2{<t))\\xi2{o) < C{p2,q2)\\4'\\xP2{Q) ■ (2.15) 

Proof. We proceed as in [12', Proposition A. 9]. We first observe that 

J(0)(x/,,X3) = / F{xh;x3,y3)dy3+ X3, ya) dys 



where 



|3:3-J/3|>1 ■J\x3-y3\<l 
Jl{(l}){Xh,X3) + J2{<P){Xh,X3) 



F{xh;x3,y3) = I Hyh^ys) —dyh, x/i G , 3:3,^3 gM. 

7r2 [xh -yhr + [xs - ysr 

For any a G M, let fa{yh) = {0? + |y/ip)^^. Then fa G L^(M^) for any r > 1 and any a 7^ 0, and 
there exists > such that 

ll/a||L'-(]R2) < „\ • 
a r 



Moreover, we have F{- ; X3, 2/3) = </>(•, ys) * fxs-ys by construction. Thus, if we take 1 < p, (7, r < 

Cr||</'(-,2/3)|||LP(R2) 



00 such that l + | = i + i, we obtain using Young's inequahty 



l^(- ;a:3,y3)||L'?(K2) < \\4'{-,y3)\\LP(R2)\\fx3-yJ 



< 



\x3 - yap" 



To estimate Ji (</>), we choose p = Pi, q = qi- In view of ()2.13p . the corresponding exponent 
r = ri satisfies 2 < ri < 00, so that 2 — G (1,2]. By Minkowski's inequality, we thus find 



ll'/i('/')(-,a;3)||i,i(K2) < / ||F(- ;x3,y3)||i9i(iR>2) dy3 < C(ri) sup ||0(-,y3)||iPi 

Taking the supremum over X3 G M, we obtain (|2.14|) . Similarly, to bound J2i<P), we take p = P2, 
q = q2. Then 1 < r2 < 2, so that 2 - ^ G (0, 1). We thus obtain 

||'/2('/')(-,a;3)||iq2{M2) < / ||i^(- ;a:3,y3)||L92(R2) dy3 < C(r2) sup ||(/)(-,y3)||iP2 

and (j2.15p follows. Finally, the uniform continuity of Ji (</>)(-, 3^3) with respect to X3 can be 
verified exactly as in the proof of [12, Proposition A. 9]. □ 

As an immediate consequence, we obtain the following useful statements. 

Proposition 2.3 Let G X{m) for some m G (l,oo]. Then J{(l)) G X'?(0) for all q G (2, 00), 
and there exists a positive constant C = C{m,q) such that 

\\Jm\xm < CUWxim) ■ (2.16) 
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Proof. If m > 1, we recall that X{m) ^ ^^(O) for all p G [1, 2]. Thus we can apply Lemma [2 .21 
with pi = 1, p2 = 2, and qi = q2 = Q ^ (2, oo), and the result follows. □ 



Corollary 2.4 Let (pi, (j)2 G X{m) for some m G (l,oo]. Then (piJ{(j)2) S XP{m) for all 
p G (1,2), and there exists a positive constant C = C{m,p) such that 

XP{m) < C\\(l)i\\x(m)\\4'2\\x{m) ■ (2-17) 

Proof. We proceed as in [I2l Corollary A. 10]. Let p G (1,2), and take q G (2,oo) such that 
I = I — ^. For any X3 G K, we have by Holder's inequality 

\\4>l{-,X3)J{4>2)(.-,Xs)\\LP(m) = ( [ Pm{\xh\'^T^^\4>l{xh,X3)\P\J{(t)2)(.Xh,X3)\P dXh) 

<( Prn{\xh\'^)\ct>l{xh,X3)\^ dXh] ( [ J(02) (a^fc, 2^3) |^ dx/^) 
= \\M-,X3)\\L2(m)\\J{<t>2){-,X3)\\Li{0) ■ 

Taking the supremum over X3 G M and using Proposition 12.31 we obtain (j2.17p . Finally, it is 
clear that the map X3 0i(', 2:3) J((/)2)(-, X3) is continuous from M into U'{m). □ 

We conclude this section with an estimate on the linear operator (jl.lip which will be needed 
in Section HI 

Lemma 2.5 Let p G [1, 2] and 2 — | < m < 00. For any (3 G N^, there exists C > such that 

||afAa;||xPM < C ||5|a;||xPM. (2.18) 

l/3|<l/3| + l 

Proof. It is sufficient to prove (j2.18p for /3 = 0. The general case easily follows if we use the 
Leibniz rule to differentiate Aw (we omit the details). 

Assume thus that oj belongs to W{m), together with its first order derivatives. Since the 
function defined in (|1.9p is smooth and bounded (together with all its derivatives), it is clear 
that 

||([/^,V)a;||xPM + ||(t^,V)^^ 

llXP(m) (m) • 

I/3|<1 

We now estimate the term {u,V)G = {K3D * w, V)G, using the fact that \K3£) * a;| < J(|w|). 
Since |w| G X^{0)riXP{0) by assumption, we can apply Lemma [2^2] with pi = 1, qi = 00, p2 = p, 
and q2 G {p, By Holder's inequality, we easily find 

II Jl(|tj|)|VG|||xP(m) < C|| Jl(|tj|)||x°°(0) < C'lll^lllxi(O) < C'll^llxP(m) , 

ll^2(|^|)|VG||| 

XP{m) < C\\J2 (1^^1)11x92(0) ^ C'|||a;| llxp(o) < C'll'^llxp(m) ■ 

We conclude that \\{u,'V)G\\xp(m) ^ II (-f^3D * V)G||xp(m) ^ C||^llxp(m)- In a similar way, com- 
muting the derivative and the convolution operator, we obtain the estimate ||(G, V)n||xp(m) < 
\\{G,'V){K3D *io)\\xp(m) < C'||Va;||xp(m,)- This completes the proof. □ 
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3 The vectorial 2D problem 



In this section we study the hnearized equation dtuj = {L — ah)uj in the particular case where the 
vorticity oj does not depend on the vertical variable. As was explained in the introduction, this 
preliminary step is an essential ingredient in the linear stability proof which will be presented 
in Section HI 

If 9^30; = 0, then £30; = 0, and the expression (|2.ip of the linear operator L becomes 
significantly simpler. On the other hand, we know from (II. lip that 

Ao; = Aio; - A2a; + K^uj - A4W , (3.1) 

where 

Aiw = {U^,V)lo = {U^,Vh)io , A3W = (n,V)G = K,V,,)G , ^^^^ 

A2UJ = {uj,V)U^ = {uh,Vh)U^ , A4W = (G,V)n = gd^.u . 

Here u = K^d * u; is the velocity field obtained from uj via the three-dimensional Biot-Savart 
law ()1.6p . Since dx^co = 0, we have dx-^u = 0, hence A^uj = in our case. Moreover, it is easy to 
verify that u = {uh, u^), where Uh = K2D *^3- Thus, we see that 

(L-aAV = := f (^'^ " D-'^ " "^^M ^ f^^^^^A , (3.3) 

^ ^ V A,W3 - a(Ai + A3)a;3 J ^ ' 

where K2UJh = {^h, ^h)Uh ^'"d A3a;3 = (K2D ^^^3, ^h)g- 

For any a € M and any m € (l,oo], the operator defined by (|3.3p is the generator of 
a strongly continuous semigroup in the space L^(m)^. This property can be established by a 
standard perturbation argument, see Lemma [3.21 b elow . Our main goal here is to obtain accurate 
decay estimates for the semigroup e*-^" as t — )• cxd. As is clear from p.3p . the evolutions for 
ujh and W3 are completely decoupled, so that we can consider the semigroups e*-^"'-'' and e*-^"'^ 
separately. The main contribution of this section is: 

Proposition 3.1 Fix m € (1, 00], a € M, /i € (0, |), and take rj E (0, ^] such that 1 + 2r/ < m. 
Then there exists C > such that 

l|e*^"'''W/i||L2(m)2 < C e~^^\\uJh\\L^{m)'2 
||e*-^"^3^3||L2(m) < Ce-'^^WusWmm) , 

for all UJ G L'^{m)'^ x Ll{m). 

Estimate (j3.5p was obtained in |1H Proposition 4.12] for m < 00, and the proof given there 
extends to the limiting case m = 00 without additional difficulty. Remark that the decay rate 
e~^^ is obtained using the fact that u;3 € Lf^{rn): If we only assume that oj^ € L'^{m) for some 
m > 1, then (j3.5p holds with rj = 0. Note, however, that uj is not assumed to be divergence- free 
in this section. 

From now on, we focus on the semigroup e'''^"'^, which has not been studied yet. To prove 
(j3.4p . we use the same arguments as in |1H Section 4.2]. We first establish a short time estimate: 

Lemma 3.2 Fix m € (l,oo], a G M, and T > 0. There exists C = C{T,m, \a\) > such that 

sup (||e*^"''^tjft||i2(„)2 +a(f)5||V,,e*^--^a;/,||i2(„)4) < C\\u;h\\L2(m)2 , (3.6) 

0<t<T^ ' 

for all ojh G L'^{rrif . Here a{t) = 1 - e~*. 



t > , 
t > , 



(3.4) 
(3.5) 
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Proof. Given € L'^{m)'^, the idea is to solve the integral equation 



tG[0,T] , 



(3.7) 



by a fixed point argument in the space Xt = {uj^ € C([0, T], L^(m)^ | < 00} defined by 

the norm 

ll^hllxT = sup \\uJh{t)\\L2(rn)2 + SUp a(t) 2 || V/^CJ/j (t) ||i2 („)4 . 
0<t<T 0<t<T 

From (|2.9p we know that ||e*^'^''~i^a;5J||xT < C'i||a;°|[ ^2(^)2, for some Ci > independent of T. 
To estimate the integral term in (j3.7p . we first observe that the velocity field U'^ defined by 
(II. 9p satisfies 



sup il + \xh\)\U'-ixh)\+ sup {l + \xh\y\VhU'^ixh)\ < 00 



In view of the definitions ()3.2p . we thus have 

11(1 + |x/i|)AiW/i|[^2(„)2 < C||V/ja;/j||i2(^)4 , 

11(1 + |x/i|)^A2Wh||L2(m)2 < C'||tJ/i||i2(„)2 . 

Using these estimates together with (j2.9p . we can bound 

ft 



(3.8) 



(3.9) 
(3.10) 



e(*-^)(^'^-2)(Ai - A2)c^/,(s)ds 

ft 



L2(m)2 



< c 







l|w/l(s)llL2(m)2 + II V/itJfe(s) 11^2(^)4 

< C||a;/,||x^ / e-i(*-^)a(s)~5ds < Ca(r)5 ||a;;,||^^ 



ds 



In a similar way, 
ft 



V 



3(*-'^)(^'>-f)(Ai - A2)t^/,(s)ds 







L2(m)4 



(3.11) 



< C 



e 2 



|(t-s 



a{t 



t-s] 

—t( l|w/i(s)llL2(m)2 + II Vfetjfe(s) 11^2(^)4 1 ds < C||w/,||xt 

S)2 ^ ^ 



Summarizing, we have shown that ||a;/j||xr ^ C'i||a;2l|L2(m)2 + C'2|ci|a(r)"'^/^||a;/i|lxr, for some 
positive constants Ci,C2- If we now take T > small enough so that C2|a|a(T)^/^ < 1/2, 
we see that the right-hand side of (|3.7p is a strict contraction in Xt- We deduce that (|3.7p 
has a unique solution, which satisfies ||ci;/i||xt < 2Ci||w°||/,2(m)2. Since uihit) = €^■^"''^00^ by 
construction, this proves ()3.6p for T sufficiently small, and the general case follows due to the 
semigroup property. This concludes the proof. □ 

We next consider the essential spectrum of the semigroup e*-^'*-'^, and begin with a few 
definitions. If j4 is a bounded linear operator on a (complex) Banach space X, we define the 
essential spectrum cTessi^ \ X) as the set of all 2; G C such that A — z is not a Fredholm operator 
with zero index, see [T7] or [Sj. The essential spectral radius of A in X is given by 



5(^;X) = sup||z| ; z G (7ess(A;X)| < 



00 . 
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If \z\ > ress(^;^), then either z is in the resolvent set of A, or z is an eigenvalue of A with 
finite multiplicity, see [5l Corollary IV. 2. 11]. In the latter case, we say that z belongs to the 
discrete spectrum of A. 

In what follows, we consider the linear operator J^a^h as acting on the complexified space 
L'^{mf', i.e. the space of all ojh : ^ ^ such that ||a;/j||j;^2(„)2 < oo. Our first result shows 
that the essential spectral radius of the operator e^^"'*^ in L'^ijnf' does not depend on a. 

Proposition 3.3 Let m S (1, oo] and q S M. Then for each t > we have 

Proof. Since ^o,h = ^h~^: the last equality in (j3.12p follows from \10\ Theorem A.l] if m < 00. 
If m = 00, then e*^*^ is a compact operator for any t > 0, hence re<js(e*-^o,h ;_L^(oo)^) = 0. To 
prove the first equality in (j3.12p . we fix t > 0. Our goal is to show that the linear operator 
Aa{t) = 6*-^°"'^ — e*''^''~2) is compact in L^(m)^. By Weyl's theorem, this will imply that both 
semigroups have the same essential spectrum, hence the same essential spectral radius. In view 
of (|3.7p we have, for all ujh € L'^{m)^, 

A«(tK = -a / e(*-^)(^'^-2)(Ai - A2)e"-^"^'^w/,ds . (3.13) 
Jo 

Let w{xfi) = 1 + \xh\- If m < 00, it follows from ()2.9p and definition (jl.l3p that 

Wwe^^f'UJhWLHm)^ < C\\e*^''^h\\L\m+l)2 < CWwujhh^m)^ , (3-14) 

for ah ojh G L^{m)^ and ah t > 0. If m = oo, we know from [13, Proposition 2.1] that 
w{—Ch + 1)"^''^ is a bounded operator in L^(oo)^, and since Ch is the generator of an analytic 
semigroup we easily obtain 

\\we^'^''OJh\\L^(mY - C\\{-^h + '^f''^e^'^''^h\\L^{ra)^ < W^hh^m)-^ , (3.15) 



for all t > 0. Now, starting from ()3.13p and using either (I3.14p or ()3.15p together with (13.9 
p.lOp . and Lemma 13.21 we find 

rt g-f(i-s) 

\\W Aa,it)uh\\L^(m)2 < C\a\ ^^^_^-^i/2 {\W''^'"-'"^h\\mm)2 + || Vfee^^^-^CJfe || ^^2 (^)4 j ds 



e-|(*-«) 



< C7|a|||a;,||^2(„)2y^ „(i_,)i/2«(,)i/2 < CIoIIKIIl^m^ • 



Moreover, proceeding as in (jS.lip . we find \\V fiA.a{t)ijJh\\ ^2 (^^^^4 < C\a\\\ujh\\L2(^^^2. Thus we have 
shown that wAa{t) and VhAa{t) are bounded operators in L'^{m). By Rellich's criterion, we 
conclude that Aa{t) is a compact operator in L^(m)^, for any t > 0. This completes the proof. 

□ 

In view of Proposition 13.31 the spectrum of the semigroup e^^°'-*^ outside the disk of radius 
e~^~^^'^^ in the complex plane is purely discrete. By the spectral mapping theorem [5], to control 
that part of the spectrum it is sufficient to locate the eigenvalues of the generator ££ci,h- Thus 
we look for nontrivial solutions of the eigenvalue problem 

^a.h^h = Aw/i , (3.16) 

where coh S L^{m)^ and A S C satisfies Re A > — y — 1. The following auxiliary result shows 
that the eigenfunctions LOh always have a Gaussian decay at infinity. 
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Proposition 3.4 Let m G (l,oo) and a € M. If ojh ^ L?'{m)'^ is a solution of ()3.16p with 
ReA > - f - 1, then ujh G L^(oo)^. 

The proof of Proposition 13.41 is postponed to Section 16.21 below. Note that a similar result 
for the nonlocal operator =Sfa,3 has been obtained in [11^ Lemma 4.5], and plays a key role 
in the derivation of estimate (jS.Sp . Thanks to Proposition 13. 4| we only need to control the 
eigenvalues of ^a,h in the Gaussian space L?'{oo)'^ . This is the last important step in the proof 
of Proposition 13.11 

Proposition 3.5 If \ is an eigenvalue of ^a,h in L^(oo)^, then ReA < — |. 

Proof. Assume that oj^ G L^(oo)^ is a nontrivial solution of the eigenvalue problem (|3.16p . for 
some a G M and some A G C Using (13. 3p . we thus have 

3 

\uh = ChUJh - -j^^h - aiUh, ^h)^h + a{^h,^h)Uh , (3.17) 

where the velocity field U'^ is defined in (jl.9p . Since ^a,h is a relatively compact perturbation 
of ^o,h = ^/i — |i both operators have the same domain, and it follows that belongs to 
the domain of Ch- In particular, we have Vhi^h ^ L^(oo)^ and \xh\ujh € ^^(oo)^, see e.g. [T3l 
Section 2]. 

In the rest of the proof, we denote by (•, •) the inner product in the complexified space 
L^(oo)^, namely 

where p{xh) = PooHxhl"^) = e'^'^l^/^. We also denote = {ujh,0Jh)- We recall that £h is a 

selfadjoint operator in L^(oo)^ which satisfies —Ch > on L2(oo)2 and -Ch > 1/2 on Lg(oo)2. 
For later use, we observe that the (unbounded) operator Uh {Uf^ ,Vh)<^h is skew-symmetric 
in L^((X))^, because the vector field p{xh)U'^{xh) is divergence-free. 

We now take the inner product of ()3.17p with ujh, and evaluate the real part of the result. 
Using the skew-symmetry of the operator {U^, V/^), we easily obtain 

ReAllw,,!!^ = {ChUJh,c^h) - -^Wo-Jhf + OiRe{{u}h,Vh)Uh,ujh) (3.18) 
= (Aw/j, w/i) - ^llw/ill^ 2aRe / p{xh){xh ■ (^h){x't ■uJK){u'^y {\xh\'^) dxh , 

where u^(r) is defined in (jl.Op . On the other hand, it follows from (I3.17P that the scalar function 
x/i • cj/i G L^(oo) satisfies 

Ax/i • w/i = Ch(xh ■ uJh) - 2xh ■ uJh - a{U^,Vh){xh ■ Wh) - 2Vh ■ (^h ■ 

Thus, proceeding as above and using the same notation (•, •) for the inner product in L^(oo), we 
find 

RQ\\\xh ■ ujh\?' = {Ch{xh ■ uJh),Xh ■ UJh) - 2\\xh ■ uJhf -2Re{Vh ■ uJh,Xh ■ ujh) ■ (3.19) 
Finally, the two-dimensional divergence • coh G Lq{oo) satisfies 

XVh-ujh = Ch{Vh ■ UJh) -Wh-uJh- a{U^, Vh){Vh ■ ujh) , (3.20) 
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hence 

ReX\\Vh-ojhf = {Ch{Vh-uJh),yh-^h) - \\^h-i^h\? ■ (3.21) 

Since V/i • w/i G Ll{oo), it follows from (|3.2ip that ReA ||Vft • ujhW^ < -|||V/i • w^Jp. Thus 
we must have ReA < — |, unless \/h • = 0. In the latter case, we deduce from (13.191) that 
ReA \\xh ■ < — 2||x/i • w/ilP, hence ReA < — 2 unless Xh ■ ^h = 0. But if this last condition 
is met, it follows from (I3.18P that ReA < — |||a;/i|Pi hence ReA < — | because oJh is not 

identically zero. Summarizing, we conclude that ReA < — | in all cases. □ 

Remark. Actually the conclusions of Proposition 13.51 can be slightly strengthened. First, in 
the invariant subspace where V/i • w/i = 0, one can show that all eigenvalues of ^a,h satisfy 
ReA < —2. This follows from the proof above if we use in addition the fact that ujh € L'q{oo)'^, 
due to the divergence- free condition. The result is clearly sharp, because if g{xh) is defined by 
(jl.Sp it is easy to verify that the function = xj^g{xh) satisfies ^a,h^h = —"^^h for any a € M. 
On the other hand, if ujh is a solution of (I3.17P such that V/i • w/i 7^ 0, we have ReA < — | if 
a 7^ 0. This follows from ()3.2ip . because we know from [101 Appendix A] that 

{Ch{Vh-^h),^h-^h) < -^\\^h-(^hf , 

unless V/i • uJh = {aixi + a2X2)g{xh) for some oi, 02 € C. But this ansatz is not compatible with 
(I3.20p if a 7^ 0. In fact, using the techniques developped in [24] or [9], it is possible to show 
that, given any M > 0, the eigenvalue equation (|3.2Up restricted to the orthogonal complement 
of the space of all radially symmetric functions in L^(oo) has no nontrivial solution such that 
ReA > — M, if |q| is sufficiently large depending on M. 

It is now easy to conclude the proof of Proposition 13.11 As was already mentioned, we only 
need to prove that estimate p.4p holds for any fi < 3/2. If pa{m) > denotes the spectral 
radius of the operator e^"-'' in L^(m)^, this is equivalent to showing that log Paim) < —3/2, see 
[5i Proposition IV. 2. 2]. But that inequality follows immediately from Propositions 13.3 1 13.41 and 
13.51 since m > 1. The proof of Proposition 13. II is now complete. □ 



4 Linear stability 

Equipped with the results of the previous section, we now study the linearized equation ()1.22p 
in its full generality. Using Proposition 12.11 and a perturbation argument, it is not difficult to 
verify that the linear operator L — qA generates a locally bounded semigroup in the space ^{m) 
for any a G M and any m E (1, 00], see Proposition 14.21 below. The goal of this section is to show 
that the semigroup e*^^"'^'^) extends to a bounded operator from W{m) to X(m) for any t > 
and any p G [1, 2], and satisfies the following uniform estimates: 

Proposition 4.1 Fix m € (l,oo], p € [1,2], a G M, and take p € (1, |), rj E (0, ^] such that 
2p < m + 1 and 2r] < m — 1. For any /3 = (/3i, /32; /^s) G N^, there exists C > such that 

||5^(e*(^--^)a;o)/.||xM2 < ^— ^||a;o||x.(„.) , (4.1) 

a{t)p 2+2 

||af(e*(^--^)a;o)3||xM < t^Tm H^ollx^M , (4.2) 

a{t)p 2+ 2 

for any ujq G W{m) and all t > 0. Moreover, V • (Jq = 0, then V • e'^^^""^^ ujq = for all t > 0. 
The proof of this important result is divided into several steps. 
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4.1 Global existence and short time estimates 

We first prove that the hnearized equation (11.22P has a unique global solution in X(m). 

Proposition 4.2 Fix m G (1, oo], p e [1, 2], and a G R. Then, for any wq € XP{m), Eq. ()1.22p 
has a unique solution u G L^^(R4.;X(m)) n C([0, oo); X^^^(m)) loii/i initial data ujq. Moreover, 
for any j3 G N^, there exist positive constants Ci, C2 (independent of a) such that 



\\dx'^{t)\\x{m) < 1 1 iffi ll^ollxp(m) , forO<t< 



a{t)p' 



a 



(4.3) 



where a{t) = 1 — e *. Finally, i/ V • cjq = 0, then V • uj{t) = for all t > 0. 



Proof. We proceed as in the proof of Lemma 13.21 Let e*^ be the semigroup generated by L, 
which is given by the explicit expression (j2.5p . The integral equation corresponding to (jl.22p is 



uj{t) = e^^ujQ-a / e(*~")-^Aw(s) ds =: {Fuj){t) , t>0 . (4.4) 
Jo 

Given A: G N \ {0} and a sufficiently small T G (0, 1], we shall solve (j4.4p in the Banach space 

Vk,T = {a;GL^,((0,r);X(m))nC([0,r];Xf„^(m)) 
equipped with the norm 



kj T < 00 



^ ( sup a(t)p 



i_i^Ml 
2^ 2 



l/3| 



l/3|<fc 



^o<i<r 



X(m) 



+ sup a(i) 2 \\d^u} 



0<t<T 



'Mi ' 



where a{t) = 1 — e~*. If (x>o G W{m), we know from Proposition 12.11 that the map t 1— t- c^^ujq 
belongs to V^^t for any T > 0, and that He^-^woUfc^T < Ci||6t;o||xp(m) foi^ some Ci > depending 
only on fc, m, p. 

Given a; G U^^t^, we now estimate the integral term in ()4.4p . Using Proposition 12.11 and 
Lemma [231 '^^ find 

||d^e^ Aw(s)||x(m) < i_i |g| < ITIIM 

a{t-s)p 2+2 2+2 

C||^llfc,T 



< 



a(t - s)p 2+ 2 0(5)2 



< s < t . 



(4.5) 



Similarly we have \\dxe^^ '^^■^Aa;(s)||xp(m) ^ Ca(t — s) 
particular case where /3 = 0, it follows that 



/9 _ !_ 

2 a(s) 2||(^||jr.y for < s < i. In the 



Auj{s) ds 



X(m) 



< Ca(t) f ||a;||fc,r , 



(*-")^Aa;(s)ds < Ca{t)^\\oj\\k,T , 



XP(m) 



< t < r . 



(4.6) 
(4.7) 



Assume now that 1 < |/3| < k. If /3' < /? and = — 1, we have from Lemma [27 

C||^lk-,r 

l/3| = l/3| 



Aa;(s)|| 



X{m) < C 2Z ll^^^('5)llx(m) < i_i ,M 

a(s)p 2+ 2 



< s < T . 
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Thus, writing Sfe^*"^)-^ = d^~^' d^' e^^-'^^ = 5f e(^4^-^3)*e(*-*)-^af', and using Proposi- 
tion [2?T] again, we obtain 



C\\io\ 



k,T 



8)2 a{s)p ^ 2 



a{t 



(4i 



for < s < t. Similarly, we have \\dxe^^ ^'^^Auj{s)\\xp(m) < a{t — s) 2 a{s) '2' Combining 
(|4.5|) and (j4.8p . we obtain the following estimate 



5f / e(*-^)^Aw(s) ds 



X(m) 



< c 



P 2 2 (i(s) 2 ds + 



a(t — s) 2 a(s) 



.1 I 



2 ds llo; 



(4.9) 



< Ca(t) p 2 ||w||fc,T 



< t < r , 

which generalizes (|4.6p . Similarly, the generalization of (14. 7p is 

3(*-^)^Aa;(s) ds 



XP(m) 



1 _]£t 

< Ca(t)2 2 



< t < T . 



(4.10) 



Summarizing, we have shown that the linear map F defined by (14. 4p satisfies the estimate 

||i^w||fc,T < C'illwollxp(m) + C'|a|T2 ||a;||fc_T , if 0<T<1, 

where C > depends only on k, m and p. Arguing as in \12\ Corollary A. 7 and Remark A. 8], 
it is also straightforward to verify that Flo € C([0, T]; X^^^(m)) if w € Vk,T- If we now assume 
that T < 6*2(1 + |ap)~^, where C2 = 1/(4(7^), we see that F is a strict contraction in l]k,T- 
As a consequence, the integral equation (|4.4p has a unique fixed point co € Vk,T, which satisfies 
||'^IU,T < 2Ci||ciJo||xp(m)- This proves that equation (jl.22p is locally well-posed in XP(m), and 
since the local existence time T is independent of the initial data, the solutions can be extended 
globally in time. Finally, since both operators L and A preserve the divergence-free condition, 
it is easy to check that, if V • a;o = 0, then the solution w of ()1.22p satisfies V ■ uj{t) = for all 
t > 0. This completes the proof. □ 



4.2 Decay estimates for the vertical derivatives 

Proposition 14.21 shows that the linearized equation (jl.22p is globally well-posed in the space 
X(m) for m > 1, but does not provide accurate estimates on the solution u}{t) = e**^^~"^^a;o 
for large times. In this section, we focus on the derivatives of uj{t) with respect to the vertical 
variable X3. Using identity (ll.23p . we shall show that d^^oj{t) decays exponentially as t ^ 00, 
provided A: G N is large enough depending on |a|. Albeit elementary, this observation plays a 
crucial role in the proof of Proposition 14.11 because it will allow us to simplify the study of the 
semigroup e*(^~"^) by disregarding most of the terms involving a vertical derivative. 



Proposition 4.3 Fix m G (1, cxd]. There exist positive constants C3, C4 such that, for all a € M, 
all A; G N, and all (Jq G X(m) with d^^uo G X(m), the following estimate holds: 

119,^3 e*(^-^)u;o||xM ^ Cse('''(\"\'+'^-'^d',.^u;oh^^) , t>0. (4.11) 
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Proof. In view of (jl.23p . it is sufficient to prove (|4.1ip for k = 0. If € X(m), we know from 
Proposition 14.21 that there exist constants Ci > 1 and C2 > 0, depending only on m, such that 
the solution u}{t) = e^^^'^^^uio of (|1.22p satisfies ||w(t)||x(m) ^ C'i||a;o||x(m) foi" t € (0,to]) where 
to = C'2/(|cKp + !)• Using the semigroup property, we can iterate this bound, and we easily 
obtain 

11 t(L-aA), , II ^ r< cC4{|Q!p+l)t II, , II 

where C3 = Ci and C4 = log(Ci). This concludes the proof. □ 



4.3 Decomposition of the linearized operator 

Motivated by Proposition 14.31 we now decompose the linear operator L — aA as follows: 

L- aA = + C3- aH , (4.12) 

where is defined in (jS.Sp and C3 = d^,^ — x^dx^ ■ We recall that the operator does not 
involve any derivative with respect to the vertical variable 2:3, and does not couple the horizontal 
and vertical components of w = (cj/j, wa)"''. In view of (|3.ip - (|3.3p . the last term in (j4.12p has the 
following expression: 

= A3 - A3 - A4 , 

where A3, A4 are defined in (|3.2p and A3 after (13. 3p . More explicitly, we have 

Hu. = (I'''] = ^~'^'''''*J''"\r, ^ (4-13) 

\H3U)J \{K3D*(^ - K2D*^^3,^)9 - 9{K3D*da:.^Uj)3j 

where K3£), K2D are the Biot-Savart kernels (11.61) . (II. 7p . and g is defined in (II. 8|) . Here * denotes 
the convolution with respect to the horizontal variables, so that 

{K2D-^^^-i){Xh,X3) = / K2D{xh - yh)^^3{yh,X3)dyh . 

Thus, unlike the operator H involves vertical derivatives, and couples the horizontal and 
vertical components of cj. As was already observed in Section [3l we have Huj = whenever 
dx^uj = 0, see Proposition 14.51 below. 

Let Ra{t) denote the semigroup generated by the linear operator + ^^3. In analogy with 
(|2.5p . we have the following representation: 

{R^{t)u){x) = -j=== j^e 2.(2t) [e'^^uj{.,y3)){xh)dy3 , t>0, (4.14) 

where a{t) = 1 — e~* and e*-^" is the semigroup generated by J^a- Since Ra{t) does not couple 
the horizontal and vertical components of a;, we can write 

where Ra^hif) and Ra,3{t) are the semigroups generated by ^a,h+^3 and ^a,3+^3, respectively. 
Using the results of Section El we obtain the following estimates: 

Proposition 4.4 Fix m € (1, 00], a G M, ^ € (1, |), and take rj G (0, |] such that 2r] < m — \. 
Then there exists C5 > such that 

(4.15) 

\\Ra,3{t)uJ3\\x{m) < C's e"''* ||a;3 |U(m) , (4.16) 
for all oj € X(m) and all t > 0. 
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Proof. Both estimates follow from the representation (|4.14p . Proposition 13.11 and estimate 
()2.1ip . The calculations are straightforward, and can be omitted here. We just remark that, 
even if V • w = 0, the map Xh i-> ujh{xh-,xz) usually has a nonzero divergence for all values 
of X3 G M. This is why Proposition 13.11 hence also Proposition 14.41 was established without 
imposing any divergence-free condition. □ 

We conclude this section with a useful bound on the linear operator H. 

Proposition 4.5 Fix m G (l,oo] and 7 G (0,1). There exists Cq > such that, for all uj G 
X(m) with dx^uJ G X(m), one has 

\\HhUj\\x{m)^ < CeWdxMkim) > (4-17) 

l|-f^3'^IU(m) < CQ{\\dx3Uj\\x{m) + \\^h\\x(^„,-j2\\dx3^^h\\x(Jn)2) ' (4-18) 

Proof. We use the expression ()4.13p of the linear operator H. Since Sxa"^ G X(m), we know 
from Proposition 12.31 that dx^u = K^d * dx^uj G X^(0). Thus, using Holder's inequality, we 
obtain 

\\gdx3u\k{m) < \\dx3u\\xH0)( [ Pm{\xh\^fg{xh)'^dXh) ^ < C\\dx3Uj\\x{m) ■ 

In particular, we have \\HhUj\\x(m)'2 ^ C'l!'9a;3'^||x{m)- 

We next consider the two-dimensional vector / = {K^d * uj — K2D *^3)h and estimate the 
term {I,'Vh)g. Using the definitions ()1.6p . (jl.7p . it is straightforward to verify that I{x) = 
Ii{x) + hix), where 

hix) = 7- / ^ ^^il {^siVh, ys) - (^siVh, X3))dy , 
4vr J^3 \x -y\-^ 

hix) = -7- / i^hiVh, y3) - ^hiVh, x^))^ dy . 

4vr J^z \x - y\^ 

Since Vhg{xh) = -g{xh)xh/2 and \xh ■ {xh - yh)^\ < \xh\\xh - yh\^'''\yh\°' for any a G [0, 1], we 
can bound 

\ih,'^h)g{x)\ < Cg{xh)\xh\ / -. r^\ujs{yh,y3) - uj3{yh,X3)\dy 

+ Cg{xh)\xh\ / -——^\ujs{yh,y3)-^^3{yh,X3)\dy. 

■J\x3-y3\<i 1^ y\ 

We now proceed like in the proof of Lemma l2.2i Integrating first with respect to the horizontal 
variable y^ G and applying Holder's inequality, we obtain 

|(/i,V/i)c/(x)| < Cg{xh)\xh\ I -. 124-^ 111 • r{^3(-,2/3) -^3(-,a;3)}|[Li(iR2)dj/3 

+ Cg{xh)\xh\ I 1 ^- r||a;3(-,y3) -a;3(-,X3)||i2m2)dy3 . 

^|^'3-?/3|<i 1^3 -y3| 

Assuming < o" < m — 1, we have the estimate ||| • |'^/||li(ir2) < C'|l/||^2(„) for any / G L^(m), 
hence 

III • r{'^3(-,y3) -'^3(-,2;3)}||l1(k2) + ||w3(-,y3) -W3(-,a;3)||i2(]g2) < C|x3 - ^3 1 ||ai.3a;3 [[^(m) • 
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We conclude that 

'^S ||x(r?i) 

'^3||x(m) ) 

which gives the bound || (/i, V)5||x(m) < C\\dxgUj3\\x(m)- 

Finally we consider the term {l2,Vh)g- Using again Holder's inequality, we obtain 

1/2(2;) I < C I - — ^—r^\u;hiyh,y3) - ^hiyh,X3)\dy 

■J\x3-ys\>i F - y\ 

+ C I ^—r2\(^h{yh,y3) - (^h{yh,X3)\dy 

< C — ^||w/i(-,y3) -a;/j(-,X3)||iim2)dy3 

+ / 1 rll'^h(-,y3) - W/i(-,X3)||i2m2)dy3 . 

^l^3-y3l<i F3 -y3| 

Since L'^{m) ^ LP{R'^) for p £ [1,2], we have \\u}hi-,y3) - ^^hi, X3)\\LP(R2y2 < 2\\ujh\\x(mf and 
\\^h{-,y3) - Wfe(-,rE3)||ip(R2)2 < \x3 -y3\\\dx3^^h\\x(m)^- I" particular, for any 7 G (0,1), 

I|wh(-,y3) - Wh(-,X3)||ip(iR2)2 < 2'''\xs-y3\'^'^\\uJh\\x^^^2\\dx3^h\\^x-{m)2 ■ 

Thus we obtain 

'""^"^ ^ C\\^h\\x(^.^-j2\\dx3^h\\x(l,)2 > 



and conclude that ||(/2, V/i)5'||x(m,) < C||a;/i||]^^^p H^j^gW/iH^^^p. This completes the proof of 
Proposition 14.51 □ 



4.4 Large time estimates 

In this section we complete the proof of Proposition 14.11 Fix m € (l,oo], a G M, and assume 
that Wo £ XP(m) for some p G [1,2]. Let uj{t) = e^^^~'^^^uJo be the solution of the linearized 
equation ()1.22p given by Proposition 14.21 Take any A; G N such that k > C4(|a|^ + 1) + 1/2, 
where C4 is as in Proposition 14.31 and choose to > small enough so that estimate ()4.3p holds 
for all t G (0,to] and all /3 G N'^ with \/3\ < k. Our goal is to control the solution u(t) for t > to 
and to establish the decay estimates (j4.ip . (|4.2p . 

To this end, we first observe that u}{t) satisfies the integral equation 

u;{t) = Ra{t - to)uj{to) - a [ Ra{t - s)Hoj{s) ds , t>to, (4.19) 

where Rait) is the semigroup defined by ()4.14p . Fix f) G (0, 1/2) such that 2f) < m — 1. By 
Proposition 14.41 we have 

ll^(t)llxM < C5e-*^(*-*o)||a;(to)||xM +C5|a| T e~^(*-^)||Fu;(5)||x(™) ds . (4.20) 

J to 

To estimate the term ||//a;(s)||x(m), we first apply Proposition 14.51 with 7 = 1/2, and then the 
classical interpolation inequality 

HQ II ^ /-(|| 11 1 — 1/A; 11 ofc ||1/^ 
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Using in addition Young's inequality, we conclude that, given any e > 0, there exists > 
such that 

Cslal ||Ftj(s)||x(™) < 4^is)\\x{m) + Ce\\d^,u}{s)\\^^^) . (4.21) 
On the other hand, since k > C4(|a|2 + 1) + 1/2, it follows from (l4lT]l that 

l|5'3^(«)llxM < C'3e-(^-*°)/2||a^3u;(to)||xM , s>to. (4.22) 
Replacing (fOT]l and K22\\ into (O0]l . we easily obtain 

Mt)hin.) < (C5llc^(to)llxM + C'MMto)\kim)) e'^^'-'°^ + ^ e'^^*"^^ 11^(5) llxM , 

for some C'^ > 0. Applying now Gronwall's lemma, and using (14. 3|) to bound ||f^(to) llx(m) ^^id 
||9^3ti;(to)||x{m) in terms of cjq, we see that \\uj{t)\\ 

x(m) ^ Ce ^ ll'^ollxp(m) t ^ to, where 
rj = f] — €. Finally, using (14. 3p again to control the solution for t < to, we conclude that there 
exists Cj > such that 

ll^(*)llx(m) = l|e*^^""^^^^o||x(m) < r::T ll^ollxp(m) ) (4.23) 

a{t)p 2 

for all t > 0. Since e > was arbitrary, estimate (j4.23p holds for any rj G (0, 1/2) such that 
2?7 < m — 1. 

To conclude the proof, it remains to find the optimal decay rates for \\ujh{t) |lx(m) ) ll'^3(^) llx(m) 
(including the value = 1/2 if m > 2), and to establish (jiJ]) . (ji^ for /? / too. First, 
combining (jl.23p . (j4.23p and using ()4.3p again for short times, we easily obtain 

Il'9:r3w(t)||x{m) = Wdx^e " -"wollxM - 1 — ll^ollxf(m) , (4-24) 

a{t)p 

for all t > 0. Moreover, if m > 2, we know from Proposition 14.41 that (|4.20p holds with f] = 1/2. 
Thus, applying Proposition 14.51 to estimate ||-fra;(s)||x(m) and using (I4.23j) . (I4.24p . we find that 
||tj(t) ||x(m,) decays like e"*/^ as t ^ oo, hence (I4.23P holds with = 1/2 if m > 2. 

Next, to obtain a faster decay estimate for the horizontal component cj/j, we use (j4.15p and 
(liTTp . Instead of (OUjl . we find 

\Mt)\\x(mP < Ce~^(*-*°)|l(a;(to))/^ltxM2 +C|a| T e-^^*"^) ||a,3^(s) ds , (4.25) 

Jto 

for any /u G (1, §)• Since ||(?2:3 6t;(t)||x(m) ^ Ce^^'^^^^^'\\dx-ji^o\\x{m) by (|1.23p . (|4.23p . we conclude 
that ||w/i(i)||x(m)2 decays like e~^* as t — )• oo, provided ;U < 1 + In other words, if /x G (1, |) 
satisfies 2/i < m + 1, we have 

\Mt)\\x{mr = \\ie'^^"'^^^o)h\\x{mr < Ge-'^*(IIH)h||x(m)2 + l|9x3^ollxM) , (4-26) 

for all t > 0. Using the arguments leading to ()4.25p and proceeding as in Proposition 14. 2^ we 
can also derive the following short time estimate, which complements (j4.3p : 

\\d^^h{t)\\x(m)2 < 1_\ \p\ (ll(^0)fellxP(m)2 + 115x3^0 llxP(m)) , < t < . (4.27) 
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Finally, to obtain decay estimates for the derivative dxOj{t), where /3 G N^, we can restrict 
ourselves to t > 2ti, where ti > is small enough so that the short time estimates (I4.3p . (I4.27P 
hold for < t < 2ti. In view of ()1.23p . we have the identity 

QP^t{L-aK)^^ = e-/^^^*-*!) 5f;>e*i(^-"^) ^{t-2t,){L-aK) ^^e^i^^-^^^WQ • 

Using the short time estimates (14. 3p . (I4.27j) with p = 2 to bound the first operator dx'l^e^'^^^~^^\ 
then the long-time estimates (|4.23p . (|4.24p or ()4.26p to treat the middle term e(*-2ti){L-aA)^ 
and finally (14. 3 p again to bound the last term dxle^^^^~'^^^ o^o, we easily obtain ()4.ip and (|4.2p . 
together with the following estimate 

l|9f(e*(^-"^)wo)h||xM2 < r^^[\\i^o)h\\xpim)^ + \\dxsU^o\kp(m)) , t > 0, (4.28) 

a{t)p 2+— 

which will also be used in the next section. This concludes the proof of Proposition 14.11 □ 



5 Nonlinear stability 

In this section we consider the nonlinear stability of the Burgers vortex and prove Theorems II. 21 
and 11.31 Our starting point is the perturbation equation (jl.lOp . which is equivalent to the 
integral equation 



2 .i 

uj{t) = e*(-^-"^)L^o + V / e(*-^)(^-"^)A^j(w(s),^j(s))ds , t > , (5.1) 



where Ni{v,w) = {K^d * v,V)w, N2{v,'w) = {v^S/jK^o * w, and is the Biot-Savart kernel 
()1.6p . We first establish the following result, which already implies Theorem 11.21 



Proposition 5.1 Fix m £ (1, oo], a € M, and take rj G (0, ^] such that 2r] < m — 1. Then there 
exist 6 = 6(0,111,7]) > and C = C{a,m,ri) > such that, for any loq £ X(m) with V ■ = 
and ||u;o||x(m) ^ Eq. (15. ip has a unique solution uj G L°°(]R+;X(m)) Ci C {[0, 00); Xiodm)), 
which satisfies 

RM^)llxM < ^^^%^e-^\ t>0, (5.2) 
ait)- 

for any multi-index /3 G N'^ of length [/3| < 1. 

Proof. Let U be the Banach space of all uj G L°°(]R+; X(m)) n C{[0, oo);Xioc{m)) such that 
V • uj{t) = for alH > and ||w||u < 00, where 



IIcjIIu = V supa(t)'2'e''iafw(t)||x(™) . 
Given ojq G X(m) such that V • wq = 0, we denote by ^> : U — )• U the nonlinear map defined by 



<i>(a;)(t) = e*(^-"^)a;o + j;i>,(^,^)(t) , t>0, (5.3) 



where <I>i, $2 are the following bilinear operators: 



^,{uj,u){t)= re(*-^)(^-"^)iV,(a;(s),cD(s))ds, i = l,2. (5.4) 
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If ||wo||x(m) is sufficiently small, we shall show that the map $ is a strict contraction in the ball 
Bk = {w G U I ||a;||u < K} for some suitable K > 0. It will follow that <I> has a unique fixed 
point CO in Bk, which by construction is the desired solution of (|5.ip . 

Since € X(m) and V -ujo = 0, Proposition 14. II shows that the map t e*(^~"^^t<;o belongs 
to U, and satisfies the estimate 



|e*(^-"^)a;o||u < Ci\\iOo\ 



(m) 



for some Ci > (depending on m, a, rj). On the other hand, if v,w £ X(m), Corollary 12.41 
implies that Ni{v,w) and N2{v,w) belong to XP{m)^ for any p £ (1,2), and satisfy the bound 

\\Nl{v,w)\\xP(m)^ + \\N2{v,w)\\xP(m)^ < C\\v\\x(m)\\Vw\\x(rn) , 

for some C > (depending on m and p). If in addition V -v = 0, then denoting u = K^u *v we 
find 

{Ni{v,v) + N2{v,v))3dxh = / Vh- {vhU3-UhV3)dxh = , (5.5) 
'■ Jk2 

for all X3 S M, hence A'^i(w,f) + N2{v,v) G W{m). As a consequence, if cj,lD € U, we have 
Nj{Lo{t),uj{t)) G X^(m)^ for j = 1,2 and all t > 0, and using Proposition 14.11 again we obtain 
the following estimate for the bilinear operators ^j: 

^y<n) ~\ Jo 



e 



2 -* p-V{t's) 



< / ^--—^\\Nj{uj{s),Cj{s))\\xP(mrds 

j=i -^0 a{t-s)p 2+ 2 

- ^ i_i M ll^('^)llx(m)l|Va;(s)||x(m) dg 

•^0 a{t-s)p 2+ 2 

< / ^ ^ds ||tj||u||a;||u < 1 ^_ ll^llull^llu • 



a(t-s)f 2+ 2 0(5)2 a(i)! 



Since we also know that Ni(u}{t),uj{t)) + N2{u}{t),uj{t)) belongs to W(m) for all t > and is 
divergence- free, we have shown that <I> maps U into U, and that there exists C2 > (depending 
on \a\, m, and r]) such that 

||^'(^)||u < C*! 11^0 lk(»n) + (^2110; llu , ||$(a;) - $(a;)||u < C2(||a;||u + HtDHu) ||w - a;||u , 

for all a;, (J € U. We now take K > such that 2C2K < 1, and assume that ||(^o||x{m) ^ ^/ (^C*!)- 
Then the estimates above show that $ is a strict contraction in the ball Bk, hence has a unique 
fixed point uj £ Bk which, of course, satisfies (jS.ip . Moreover ||w||u < 2Ci||u;o||x(m)) hence ()5.2p 
holds with C = 2Ci. This concludes the proof. □ 

Remark. The size 5 of the local basin of attraction of the Burgers vortex aG in X(m) depends 
a priori on a, m, and 77. However, as announced in Theorem 11.31 the dependence on the decay 
rate 77 can easily be removed by the following (standard) argument. Given m > 1, we first 
choose 7] = fj{m) = min(^, ^^^) and apply Proposition 15.11 with that value of rj. We thus obtain 
a constant ^ > depending only on a and m such that, for any uq € X(m) with V • cjq = 
and ||a;o||x(m) < Eq- (15. ip has a unique solution uj G L°°(R+;X(m)) n C {[0, 00); Xioc{rn)), 
which converges exponentially to zero as i — )■ 00. In particular, given any t] G (0, ^] such 
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that 2r] < m — 1, there exists T = T{r]) > such that ||w(t)||x(m) ^ ^ for t > T, where 
6 = 6{a, m, rj) is the constant given by Proposition 15.11 By uniqueness of the solution, we 
conclude that oj satisfies (j5.2p for any admissible value of r]. 

In view of Proposition 15.11 and the remark that follows, the proof of Theorem 11.31 will be 
complete once we have established the improved decay estimate (jl.20p for the horizontal com- 
ponent ojh. A convenient way to do so is to repeat the proof of Proposition 15. II using a different 
function space, which incorporates a faster decay rate as t — )• oo. Given /i G (1, 1 + 77), where 
r/ S (0, \] is as in Proposition 15.11 we introduce the space V C U defined by the norm 

A^=0,l|/3|<1 *>° 

As in the remark above, we can assume here (without loss of generality) that ||9xWo||x(m) is 
finite and arbitrarily small, for all /3 € with \j3\ < 1. Using Proposition 14.11 we thus obtain 

||e*(^-"^)a.o||v < C3 E ll^f^ollxM , 

I/3|<1 

for some C3 > 0. On the other hand, ii v,w G X(m), the following estimates hold for any 
i^G (1,2): 

\Wl,h{v,w)\\xv{m)^ < C\\v\\x(^rn)\\^Wh\\x{m)^ , 

"^3 llx(m) ll^a;3^IIx(m) ) ; 
\\dx3Nj{v,w)\\xP(m)3 < C{\\dx3v\\x(rn)\NMk{m) + \Hx{rn)\\dx3'^w\\x{m)) ■ 

We now estimate the bilinear operators ^j{uj,u}) for uj,uj S V. First, using (14.28^ . we find 
for t > 1: 

liaf <I>i,,(..,a.)(t)b(„)2 < r {e(*-^)(^-"^)iVi(a;(s),c:;(s))},||^(^)2 ds 

Jo 

rt ^-n{t-s) 

- ^ 1,1 \ii (.\\^l,h{^is),'^{s))\\xP{mp + \\dx3Ni{uj{s),Cj{s))\\xP(m)3)ds 

•^0 a{t-s)p 2+ 2 

f-t ^-fi{t-s) 

- ^ T7T-iir(ll^(^)llx(m)l|Va:'fe(s)i|(x{m))2 

•^0 a{t-s)p 2+ 2 

+ l|5x3w(s)||x(m)||Va}(s)||x(m) + ||^(s) ||x(m) Ilea's VlD(s) ||x(m)) ds 
i-t ^-fj,{t-s)^-{^i+r])s 

< C / — - — ds ||a;||v||w||v < Ce^^^llcjUvllt^llv • (5.6) 

In the last inequality, we have used the definition of the norm in V and the fact that ^+?7 < l+2r]. 
The bound ()5.6p also holds for t < 1, and can easily be established using ()4.ip instead of (I4.28p . 

Next, to bound a^3$i,/i(tj, a)), we recall that d^.^e*^^'"'^^ = e-*e*(^-°^)(9a;3- Applying (l4l^ . 
we find 

\\dxM^^,h{^,Cj){t)\\x(u.)^ < /"*e-(*-^)||af{e(*-^)(^-^)a,37Vi(a;(s),cD(s))},||^(„)2d5 







•^0 a{t-s)p 2+ 2 

/■* g-(/^+l)(t-«)g-(M+»?)« (Jg^-(P-+V)i 

- ^ / T^* ll'^lkll'^llv < ll^llvil^l 

•^0 a{t-s)v 2+ 2 a{s)2 a{t)p^ 2 
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Similarly, for A; = 0, 1, we can estimate 9^ $2,/i(^>^) as follows: 



X3 

t 





^ ^ / ^-^\\di^N2{u;is),^is))\\xr:im)^ds 

■^0 a{t-s)p 2+ 2 

/■* g-(M+fc)(t-s)g-(A'+»?)s (JQ-{^J-+kr|)t 

- ^ / TTTTM T^^ ll'^llvll'illv < i^M^fe_3 ll^llvll^llv • 

■^0 a(t-s)p 2+ 2 a(s)2 a(t)p 2 ^2 2 

Finally, using (j4.2p . we obtain for the vertical components of ^j{uj,Gj): 



■^0 a(t-s)p a"*" 2 

- ^ / TTTTM ll'^lkll'^llv < 1 ,M, fc_3 ll^llvli^llv • 

■^0 a(t-s)p 2^ 2 0(5)2 a(t)p 2 "^2 2 

Summarizing, we have shown that $ defined by (|5.3p maps V into V and satisfies the following 
bounds: 

II^MIIV <C^Y1 \\d^x^^\Wm)+CA\W\\N > 
I/3|<1 

||$(a;) - ^>(tD)||v < C4(||a;||v + ||'^||v)||w - cu||v , 

for all G V. If K = 2C3 X]|/3|<i ||c?^^o||x(m) is sufficiently small, it follows that $ is a 
strict contraction in the ball Bx = {(^ G V [ ||a;|jv < K}, hence has a unique fixed point there. 
Denoting by ijj{t) the solution of (jS.ip given by Proposition 15.11 this implies that t 1— )• uj{t + T) 
belongs to if T > is sufficiently large. In particular, io{t) satisfies (jl.20p for some suitable 
C > 0. The proof of Theorem 11.31 is now complete. □ 



6 Appendix 

6.1 Proof of Lemma 11.11 

Let X £ C^(]R^) be a cut-off function such that x{^h) = 1 if \xh\ < 1 and xi^h) = if \xh\ > 2. 
Given i? > 0, we denote xni^h) = xi^h/R), so that \VhXR{xh)\ < C / R. For any X3 £ M, we 
define 



/(^3) 



Since 0)3 E for some m > 1, it is easy to verify that ||/ — 

the other hand, for any test function ijj G Co°(M), we have 



0j3{Xh,X3)xR{Xh)dXh ■ 

I — )• as 



[ fix3)^{x3)dx, < /■/^(X3);^(x3)dx3 

Jr dx3 Jr dx3 



dV' 



dx3 



00. On 



(6.1) 



The last term in the right-hand side converges to zero as i? — )■ 00. To treat the other term, we 
observe that 



/ /ij(2;3)T^(2;3)dX3 = / 0J3{Xh,X3)xRixh)^{x3) dX3 

Jr dx3 dx3 



{^3 , 



dx3 
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where 4>R{xh,X3) = XR{xh)ipix3) and (•, •) denotes the duahty pairing of X''(M^) and Cq° 
Now, since V • a) = in the sense of distributions, we have 

^'^^ ' ^ ^ ~^^3^ ' ^^'^ ^ ^^'^ ' '^'^ ' ^^"^ ^ "('^/^ ' ^'«'^^?) ' 

SO that 

/ fRix3):\ ix3)d.X3 = - U}hiXh,X3) ■VhXRixh)ip{x3)dxhdx3 . 

Jr dxs J^3 

Using the inclusion L'^{m) ^ L^(M?) and the definition (jl.l5p of the space X{m), we thus find 



/ fRix3)P^ix3)dx3 < ^ ||tD/i||x(m) 

Jr dx3 R ' 



R-^oo 



. 



Returning to (jG.ip . we conclude that the left-hand side vanishes for all ip € C^(]R), hence 
= in the sense of distributions. Since / G BC{M.), it follows that / is identically constant, 
which is the desired result. □ 

Remark. If uj[x,t) is any solution of p.lOp that is integrable with respect to the horizontal 
variables, we can define 



(t){x3,t) = I U}3{xh,X3,t)dXh , X3 € M , t>0. 

As was observed in [12j, this quantity satisfies a remarkably simple equation 

dt(t){x3,t) + X3dx.^(j){x3,t) = dl.^(l){x3,t) , (6.2) 

which can be solved explicitly. However, if uj{-,t) G X(m)^ for some m > 1 with V • uj{-,t) = 0, 
Lemma [TT] shows that </>(x3, t) does not depend on X3, and ()6.2p then implies that (j){x3, t) is also 
independent oft. Thus, as was already mentioned, we can restrict ourselves to the particular case 
where (j) = without loss of generality. Being unaware of this simple observation, the authors of 
[12] have stated their stability result in a seemingly more general form, allowing (apparently) for 
nontrivial functions (j)(x3,t), but thanks to Lemma [1.11 (which also holds in the slightly different 
functional setting of [12] ) the simpler presentation adopted here in Theorem 11.21 is exactly as 
general. 

6.2 Proof of Proposition 13.41 

This final section is devoted to the proof of Proposition 13. 4[ which shows that eigenfunctions 
of ^a,h corresponding to eigenvalues outside the essential spectrum have a Gaussian decay at 
infinity. For the nonlocal operator ^a,3, the same result was established in [111 Lemma 4.5] 
using ODE techniques, but we prefer using here a more fiexible method based on weighted 
estimates. In fact, we shall consider a more general elliptic problem of the form 

-Cf + F{x,f,Vf) + Xf = h, xGM", (6.3) 

where the unknown is the vector- valued function / = (/i, . . . , /at)^. Here and below we denote 
by£ = A + |- V + ^ the analog of operator (12. 2p in dimension n. The data of the problem are 
the functions F : x x C"^ and /i : M" ^ C^, and the complex number A. 

For m € [0, 00], we denote by Lp'im), H^{m) the following complex Hilbert spaces on 

L\m) = {/ G l2(M", C) / |/(x)|Vm(l^l') dx < ooj , 

H\m) = {/ G L\m) | 8^ J G L^m) for j = 1, . . . , n} , 
where pm is the weight function defined by ()1.12p . Our main result is: 



pn. 
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Proposition 6.1 Let m € [0,oo), A € C, /i € L^(oo)^, and assume that F is a continuous 
function satisfying 

\F{x,p,Q)\ < A{x)\p\ + B{x)\Q\ , for all {x,p,Q) eM."" xC^ xC"^ , (6.4) 
where A and B are hounded, nonnegative functions such that 

lim sup A{x) = lim sup B{x) = . (6-5) 

\x\>R \x\>R 

//ReA > f - f , then any solution f € H^{m)^ of ([O]) satisfies f G H^{oo)^. 

Proof. The proof is a simple modification of Proposition 12], which in turn is inspired by 
a recent work of Fukuizumi and Ozawa [6] where decay estimates are obtained for solutions of 
the Haraux-Weissler equation. For k > 1, e > 0, and 6 G [0, m], we define the weight functions 

CkAx) = e^^^^^ , Ce{x) = (1 + \x\y , X G . (6.6) 



Multiplying both sides of (|6.3p by C,eik,ef and integrating by parts the real part of the resulting 
expression, we obtain the identity 

/ Ce6,.|V/|2dx + Re / /•(V(CeeM),V)/dx+ / • V(Cea,,) dx 

JR" JR" JM" 4 

= -Re [ Ce^kJ-F{xJ{x),Vf{x))dx+('^-ReX) [ Ce^kMl^ dx (6.7) 
+ Re/" Oefc,J-/idx. 



Clearly, 

^^'''^"^^ = (4fc + |xp)2 ^^'-^^) ' = TTW^'^^^ ■ ^^'^^ 

Thus, the second term in the left-hand side of (|6.7p can be written in the following way: 



Re / / • (6,eVCe, V)/ dx + Re [ / • (CeV^fc,,, V)/ dx 

= -/ |/PV.(^^^) dx + Re / /•(C.Ve,,„V)/dx 
Jr" V1 + |x|2/ 



To bound this quantity from below, we observe that 



[ \f\'^k,eX-V-^dx<29' [ ,i^|/r^dx. 
Moreover, for each r/i > 0, 
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Thus, using the expression (j6.8p of V$,k,€, we find 



(1-6)^ f k\eCk,e\^f\^ 



where C = nO + O"^ does not depend on k and e. We next consider the third term in the left-hand 
side of (I6.7p . which satisfies 

To estimate the right-hand side of (|6.7p . we use (j6.4p and obtain, for each r]2 > 0, 
-Re/ Ce^kJ-F{xJ{x),Vf{x))dx < [ Ceefe,.^l/l'dx+ / Cee^^l/llV/l dx 

< /" Ce^kJA + ^)\ffdx + rj2[ Ce4,.|V/pdx . (6.11) 
Finahy, for each 773 > 0, we have 

Re/" CeCkJ-hdx < 7J3 [ C06^,|/|2dx + -^ / CeCfc.el/il' dx . (6.12) 

JR'i JR" 4?73 jRn 

Substituting (j6.9p - (j6.12p into (16. 7p . we arrive at our basic inequahty: 

/ Ce&,. V/ dx+ / (— I dx 

Jr" Jr" (4/c+[x|^)^ V l-r]i l + |x|^/ 

^ / CeCkJ(-r^ + j-ReX + A + ^ + rjs-l)\ff + ^\h\'}dx. (6.13) 
JiRn lVl + |x|^ 4 4r/2 2/ 4r/3 J 

To exploit (j6.13p . we first take r]i = r]2 = \ and 6 = m. Using (j6.5p and the assumption 
that Re A > ^ — y , we see that there exists ii > independent of A; > 1 such that, if 7/3 > is 
sufficiently small, the following inequality holds: 

m11''m2^^2 dx < C / oa,.|/Pdx+-L /■ C.a,.|/i|'dx, 
Jr" (4A; + |x|2)^ J\x\<r Jm.^ 

where the constant C > is independent of A; > 1. Thus, taking the limit — )• 00 and using 
Fatou's lemma, we obtain 



6(1-6) 
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/ (l + |x|2)'"e^l^l'|x/pdx < C(i?) / \ffdx + ^[ (l + |x|2)'^e^l^l>pdx 
Jr" J\x\<r 4r/3 J^n 



which shows that e~8~l^l / g L2(r2) for any e > 0. Next we choose ryi = 5, 7?2 = |, % = 1, and 
= in (|6.13p . Taking again the limit A; — t- 00 and using Lebesgue's dominated convergence 
theorem, we find 

\[ e^N'lV/l^dx + i^ / e^N^|x/|2dx<C / e^l^l Vl' dx + ^ / eVN^|/i|2dx, 

" JR" 24 Jjgn Jgn 4 J]g„ 
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where the constant C > does not depend on e > 0. This inequahty shows that 

- / e^l^lV/|2dx+^ / e^N'|x/Pdx <c[ e^l^l' dx+^ / e^l^l'j/ipdx, 
8 Jr" 48 J\x\<b' 4 J^n 

for some R' > independent of e > 0. Taking now the hmit e — >■ 0, we conclude that / G H^{oo), 
which is the desired result. □ 



Proof of Proposition 13.41 We consider the eigenvalue equation (j3.16p . which can be written 
in the form 

- Ch^^h + aAiUh - aMujh + (^X + ^^Uh = , (6.14) 

where Ch is given by (12. 2p and the operators Ai, A2 are defined at the beginning of Section [3l 
We recah that |Aic<j/i| < jV/iW/jl and IA2CJ/1I < \V hUf^\\ujh\, where the velocity profile 
satisfies (|3.8p . Assume that ReA > — y — 1 and let coh G H^{m)'^ be a solution to (j6.14p . 
Applying Proposition 16.11 with n = N = 2, F(x, fjV f) = aAif — aA2/, and h = 0, we obtain 
LOh G -ff^(oo)^. This completes the proof of Proposition 13. 4[ □ 
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